Growth of epitaxial nanowires by controlled coarsening of strained islands 
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We show that elongated nanowires can be grown on crystal surfaces by allowing large strained 
two-dimensional islands to desorb by varying the adatom supersaturation or chemical potential. 
The width of the wires formed in this process is determined by a competition between the repulsive 
elastic interactions of the long edges of the wires and the thermodynamic driving force which tends 
to decrease the distance between these edges. The proposed mechanism allows for control of the 
wire sizes by changing the growth conditions, in particular, the vapor pressure of the material that 
is being deposited. 
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Strained heteroepitaxial growth of semiconductor and 
metal films has been of considerable interest in the re- 
cent years as it provides a versatile route to fabrication of 
nanostructures that have potential applications as elec- 
tronic, optoelectronic and memory devices. Well stud- 
ied examples of such structures include 3D islands or 
quantum dots in strained group IV and III-V systems. 
Mismatch strain can also give rise to elongated wire-like 
shapes with heights in the single digit nanometer range, 
since the increased perimeter of these structures leads 
to efficient relaxation of strain Q, 0, 0, 0, S ■ While 
nanowires may have application as non-lithographically 
fabricated interconnects, they can also be used to gain in- 
sight into fundamental physics of low-dimensional quan- 
tum systems. 

Although a competition between elastic and boundary 
energies of 2D islands can give rise to elongated shapes, it 
is difficult to control the island widths and lengths in an 
independent manner, since the aspect ratios are closely 
related to their area Jj. Recent experiments have also 
shown that strain in 2D islands can be also efficiently re- 
laxed by ramification of the edges Q or by formation of 
curved boundaries 0]. In this article, we show that elon- 
gated nanowires can be formed during the dissociation of 
large islands as a result of the competition between the 
repulsive clastic interactions between the edges of the 
island and thermodynamic driving force which tends to 
decrease the distance between the edges. Key advantages 
of the proposed mechanism are that 1) it allows the con- 
trol of the width of the wire by changing the chemical 
potential or vapor pressure of the adatom reservoir and 
2) elongated nanowires can also be formed in systems 
where ramified or curved boundaries are responsible for 
stress relaxation of large islands 0, ■ 

If island growth is limited by the attachment of atoms 
to its perimeter, the outward normal velocity can be writ- 
ten as 

v n (s) = k [A/j, — k((3(9) + (3" (9)) - rV • u(a)] , (1) 

where s denotes the arc-length (refer to Fig. 1), k is re- 
lated to the mobility of the perimeter, A/i gives the driv- 
ing force on a straight edge, K is the curvature (positive 
if the center of curvature is inside the island), (3(9) is 
the energy required to create a unit length of the island 




FIG. 1: Schematic of the strained island or surface domain, 
whose boundary is represented by S, the Cartesian coordinate 
and the unit-normal of a point on the boundary is denoted 
by x(s) and n(s), respectively, where s is the arc-length. In 
the case of a strained island, r represents the product of the 
equibiaxial mismatch stress and the height of the island. 

boundary, and the last term is the contribution due to 
elastic relaxation of material close to the island perime- 
ter. Here, r is the product of the equibiaxial mismatch 
stress 9] and the height of the strained island, or the 
difference between the surface stresses of the domain and 
the substrate in the case of strained surface domains 0]. 
The divergence of the displacement field can be written 
using elastic Green's functions ^3] as 

V u(s) = (1 v2)T f [*(s)-x(s')]-n(s') , ( ) 
Js |x(s)-x( S ')| 3 ' 

where v and E are the Poisson's ratio and the Young's 
modulus of the substrate, respectively, and the Cartesian 
coordinate x and the outward unit-normal n are shown 
in Fig. 1. In what follows, we study shape evolution 
of islands by discretization of their perimeter into small 
linear segments and by moving these segments according 
to the velocity in Eq. (1). The integrand in Eq. (2), 
however, is singular due to the divergent nature of the 
elastic fields at the island perimeter. It can be regularized 
using an elastic cut-off length S, which we have set to 2 
A in the numerical calculations. 

We first consider the growth of islands in a system 
with a 2-fold symmetric boundary energy, (3(9) — (3q(1 — 
e cos 29), with the anisotropy factor e = 0.1. The effect of 
strain on the island shapes can be conveniently discussed 
in terms of the dimensionless parameter a — Fj (3$, where 
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FIG. 2: (a) Growth shapes of a strained island with a 2- 
fold symmetric boundary energy. The red curve in the inset 
graph shows the evolution of the island aspect ratio (ratio 
of the lengths of major and the minor axes), while the black 
dashed curve denotes the equilibrium aspect ratio when strain 
effects are not included, (b) Dissociation shapes of the island 
marked 5 in (a). When A/j, < 0, the island initially shrinks 
in the horizontal direction before it begins to steadily grow 
along the vertical direction. All the lengths are given in nm. 

F = (1 — i/ 2 )t 2 /(ttE). In the case of surface domains 
studied in [8( and [Tj, a m 0.3, while it is smaller 
(a « 0.05) in the case of 2-fold symmetric Si(001) is- 
lands considered in §. The growth shapes (Afi > 0) of 
a circular nucleus for a — 0.2 are shown in Fig. 2(a). We 
find that the aspect ratio continously increases as the is- 
land grows in size, far exceeding the aspect ratio of the 
equilibrium shape without strain effects (refer to the in- 
set in Fig. 1(a)). The evolution of the aspect ratios is 
qualitatively similar to the experimental observations of 
island growth on Si(001) by Zeilasek et al @. We can 
understand the fast outward growth of the regions with 
smaller boundary stiffness by comparing the contribu- 
tions to the normal velocity from elastic and line-tension 
effects. While the contribution from the line energy is 
proportional to the curvature k, the elastic contribution 
is opposite in sign and can be shown to be proportional 



to Klog(K<5). Therefore, regions with large curvature and 
small boundary stiffness tend to move outwards at a more 
rapid rate compared to the regions of smaller curvature 
due to larger contribution from the elastic relaxation ef- 
fects in the former case. 

Next, we consider the evolution of the shape marked 
5 in Fig. 2(a), when the driving force A[i < 0. In this 
case, as the material close to the minor (shorter) axis 
of the island begins to desorb, the width of the island 
decreases, while the regions near the major axis show 
a very small outward movement (refer to shapes 1-3 in 
Fig. 2(b)). While the contributions of the line energy 
and elastic relaxation to the normal velocity are smaller 
than the contribution from the driving force near the mi- 
nor axis, its effect is offset to a large extent by the elastic 
contribution near the major axis due to the higher curva- 
ture of this region. As the island shape evolves, the width 
of the island steadily decreases and eventually reaches a 
point where the repulsive elastic interaction between its 
nearly straight edges (shape 4 in Fig. 2(b)) becomes com- 
parable to the driving force, A^i. At this point, the com- 
petition between these two effects leads to an optimum 
width of the island, while the top and bottom ends con- 
tinue steadily move in an outward direction giving the 
elongated shapes shown in Fig. 2(b). An animation of 
the wire formation process can be found at the web-site 
given in Ref. |12|. 

A quantitative expression for the width w of the wire 
can be derived by considering the equilibrium configu- 
ration of an infinitely long wire for Afj, < 0. Since the 
velocity of the straight edges of the wire must vanish in 
equilibrium, we have, 

/ 2F\ 2F 

v n = k A[i H 1=0, or w = — - — . (3) 

\ w J —Afj, 

In the calculations in Fig. 2(b), we have chosen A/x = 
-0.004/V<5 and F = 0.2/3 , which gives w = 100(5, in 
agreement with the observed width of the wire. Since the 
width of the wire depends on Afi, it can be controlled by 
varying the supersaturation or the vapor pressure of the 
material that is being grown. The competition that leads 
to the wire-like shapes in Fig. 2(b) also occurs in the case 
of elongated stripes of 7x7 domains of Si(lll) that coex- 
ist with the lxl domains |llj| at temperatures where the 
free energy of the former structure is larger than the lat- 
ter (which is equivalent to A^l < 0). As in Eq. J3J), the 
stripes are stabilized by the repulsive elastic interactions 
between the edges. It should also be noted that since 
Eq. ||3J) does not depend on the form of boundary en- 
ergy, elongated shapes could also be achieved in systems 
with other types of symmetry, as illustrated in the next 
paragraph. 

In Fig. 3(a), we show the growth shapes of an is- 
land with a 4- fold symmetric boundary energy, (3(9) — 
/3 (1 - 0.033 00846*), for F/f3 = 0.2. Here, as the cir- 
cular nucleus grows in size, its shape initially evolves 
to a square, after which the boundaries acquire con- 
cave features with curvatures that continuously increase 
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FIG. 3: (a) Growth shapes of a strained island with a 4- 
fold symmetric boundary energy. Note that as the island 
grows in size, its edges become more concave in curvature, 
(b) Dissociation shapes of the island marked 6 in (a). When 
A/i < 0, the island initially shrinks in size before it begins to 
steadily grow along the vertical direction. All the lengths are 
given in nm. 

with size. Such concave shapes, which lead to relax- 
ation of mismatch strain, have been very recently ob- 
served on strained domains on Si(lll) @. The dissoci- 
ation shapes of the island marked 6 in Fig. 3(a), with 



A/i = -0.008/V<5, are shown in Fig. 3(b). The shape 
evolution involves two stages: First, the island shrinks in 
size until it reaches the shape marked 3 in Fig. 3(b) at 
which point the system slowly elongates in the vertical di- 
rection (refer to the animation in Ref. [12|) due to the in- 
stability discussed in Refs. pj and 0. As the island elon- 
gates at fixed A/x, the curvatures close to top and bottom 
edges become large, leading a steady outward propaga- 
tion, while the long edges become straighter and achieve 
the optimum width given by Eq. It should be pointed 
out that the shape marked 3 is a shallow metastable state 
and the direction of elongation (horizontal or vertical) 
depends on the small noise in the simulations due to fi- 
nite arithmetic precision. In an experimental situation, 
thermal noise will select the direction of propagation. 

While most experiments to date focus on growth or 
equilibrium shapes of islands, the shapes of strained 
Ca/CaF2 islands observed by Batzill and Snowdon can 
perhaps be understood on the basis of the present work. 
While island shapes on atomically smooth CaF2 surfaces 
are consistent with the equilibrium theory in Ref. pj , the 
widths of the long wire-like shapes on stepped surfaces 
were found to be significantly larger than the predictions 
of this model. A possible explanation for the discrep- 
ancy is that the step-edges can provide favorable sites 
for attachment of Ca atoms, and therefore give rise to 
an effective A/i < 0, which would lead to coarsening of 
large islands. As the large islands lose Ca atoms, the 
competition between elastic interactions and the coars- 
ening effect can result in the observed wire-like shapes. 
Of course, more experimental and theoretical work is re- 
quired for a quantitative verification of this argument. 

In conclusion, we have developed a numerical method 
which gives growth shapes of strained domains that are 
similar to the shapes observed in recent experiments. Us- 
ing this method we have shown that straight and elon- 
gated nanowires can be grown by controlled desorption 
of large islands. 
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